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Formulas for the Percentage Points of the Distribution 

of the Arithmetic Mean in Random Samples from 

Certain Symmetrical Universes * 

By Uttam Chand 

Using the method of Fisher and Cornish, the 100e% point of the distribution of the 
arithmetic mean in random samples of size N from any universe having finite cumulants, 
K 1} K 2 , K3, . . . , is expressed to order 1/iV 2 as a function of (i), the first five positive integral 
powers of the 100e% point of a standardized normal variable, and (it) the quantities K\, 
(K 2 ) V2 , and K r /(K 2 ) r / 2 for r=3, 4, 5, and 6. The numerical coefficients involved are evaluated 
for the case of sampling from the normal, rectangular, double-exponential, sech, and sech 2 
distributions, and the accuracy of the resulting formulas illustrated by numerical examples. 



Consider y x . . . y N a random sample of size N 
from a population having finite cumulants K r = 
K r (y), (r=l,2, . . .). Let a r =a T (y)=K T /K 2 r < 2 be 
the relative cumulants, (r=l,2, . . .), and let 
X and a specify the location and scale parameters 
of the population distribution, respectively. 

Let y denote the sample mean and consider the 
standardized variate 

,_y-gi(y) 



Applying the formulas of Fisher and Cornish 
[1, 2] 2 and utilizing the relations 

K x {y)=KM 

K 2 (y)=K 2 (y)/N 

ar(0=0 



and 



a 2 (t) = l 



ar(t)=a r (y)=a r (y)/N* ,(r=3, 4, . . .), 



the one-tail e-probability point, t € , of the standard- 
ized variate t can be expressed to order of 1/N 2 as 

1 Revision of a paper written during the summer of 1947 when the author 
was a guest worker at the National Bureau of Standards, and which was read 
by title at the Eleventh Summer Meeting of the Institute of Mathematical 
Statistics, held in Madison, Wis., Sept. 7 to 10, 1948. 

2 Figures in brackets indicate the literature references at the end of this 
paper. 
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where x=x ( is a normal deviate exceeded with 
probability e. 3 

In particular for a symmetrical universe we have 

dr(y)=0 for r=3, 5, 7, . . ., 

whence 

U=X+ MN ( - r3_3a;)+ 2V 2 [72 6 (*"- 10^+16*) - 



f^ (3* 5 -24a; 3 +29. 



a?)]- 



3 Miriam L. Yevick has drawn to my attention the fact that comparison 
of the above expression for t* with the corresponding expression given by 
Simaika [3], which includes terms involving 7i(=as) and 72(=gu) only, reveals 
two errors in his numerical coefficients: the coefficient of 7! given as "—29/ 
381" should read "-29/384", and the coefficient of 7172 given as "-107/288" 
should read "+107/288". 
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It is interesting to note that for the one-tail 
0.04163-probability point, for which x 2 e =3, we 
have 

t ( =x € + 0(N- 2 ) 
as N— > °° . 

The density functions and the requisite cumu- 
lants of five symmetrical distributions considered 
explicitly are listed in table 1, and the values of 
tt computed for iV=10 from 1-term (normal 
approximation), 2-terms, and 3-terms of the 
immediately preceding formula are compared with 
the corresponding true values in table 2. 



Table 1. Five symmetrical distributions and some of their 
cumulants 



Distribution 


Density function 


K 2 

(variance) 


a 4 (=T 2 ) 


H 




-(*/-X)2 








Normal . 


1 2o-2 


0-2 










Rectangular 


1 1 ^ 2/-X ^ 1 

- for — - < < - 

cr 2 a ~ 2 


0-2 

12 


6 
5 


48 

7 






-\y-\\ 








Double exponential- . 


1 o- 

2 


2o-2 


3 


30 


Sech 


— sech ( ) 

7T0- \ a J 


7r2(r2 


2 


16 


Sech 2 


— sech 2 ( J 

2<r \ a ) 


TT 2 <T 2 


6 
5 


• 



Table 2. Comparison of various approximate and exact values of probability points, t e , of the standardized variate t 

N=10 



Distribution a 


One-tail probability levels 


0.001 


0.005 


0.010 


0.025 


0.050 


0.100 


0.250 


Normal: 
1 term (exact) . 


3. 0902 

2. 9890 
2. 9837 

2. 9825 

3. 3432 
3. 3038 
3.3114 

3. 2589 
3. 2426 
3. 2449 

3. 1914 
3. 1840 
3. 1646 


2. 5758 

2. 5290 
2. 5270 
2. 5269 

2. 6929 
2. 6839 
2. 6839 

2. 6539 
2. 6492 
2. 6494 

2. 6226 
2. 6219 
2. 6207 


2. 3263 

2. 2983 
2. 2973 
2. 2971 

2. 3965- 
2. 3941 
2. 3932 

2. 3731 
2. 3712 
2.3711 

2. 3544 
2. 3548 
2. 3540 


1. 9600 

1.9517 
1. 9515- 
1. 9515 

1. 9806 
1. 9822 
1. 9814 

1. 9737 
1. 9738 
1. 9736 

1. 9682 
1. 9691 
1. 9681 


1.6449 

1. 6473 
1. 6474 
1. 6474 

1. 6388 
1. 6409 
1. 6406 

1. 6408 
1. 6415- 
1.6414 

1. 6424 
1. 6430 
1. 6426 


1.2816 

1. 2903 
1. 2905+ 
1. 2905 

1. 2598 
1. 2612 
1. 2612 

1.2671 
1. 2677 
1. 2677 

1. 2729 
1. 2730 
1. 2732 


0. 6745- 

.6831 
.6833 
.6833 

.6530 
. 6531 
. 6531 

.6602 
. 6605- 
. 6605- 

.6659 
.6657 
.6661 


Rectangular: 
2 terms 


3 terms - - 


exact 

Double exponential: 

2 terms 

3 terms 


exact 

Sech: 
2 terms . . _ _ 


3 terms 


exact 

Sech2 : 
2 terms.. - _. __ 


3 terms 


exact 





a The exact values of t ( , for the normal distribution were taken from Kelley [4]; for the rectangular and double-exponential distributions, from an unpublished 
table prepared in connection with [5]. The exact values of U for the sech and sech 2 distributions were computed under the supervision of Irene A. Stegun of the 
National Bureau of Standards Computation Laboratory from formulas prepared by Julius Lieblein of the National Bureau of Standards Statistical Engineering 
Laboratory. 



The author expresses appreciation to Churchill 
Eisenhart for suggesting this problem and to Lola 
S. Deming and Celia S. Martin of the Statistical 
Engineering Laboratory, National Bureau of 
Standards, for the numerical computations. 
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